Theory of vortex matter transformations in the high T c superconductor YBCO 
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^■y Flux line lattice in type II superconductors undergoes a transition into a "disordered" phase 

like vortex liquid or vortex glass, due to thermal fluctuations and random quenched disorder. We 
quantitatively describe the competition between the thermal fluctuations and the disorder using the 
Ginzburg - Landau approach. The following T — H phase diagram of YBCO emerges. There are 
O^l ' just two distinct thermodynamical phases, the homogeneous and the crystalline one, separated by a 

single first order transitions line. The line however makes a wiggle near the experimentally claimed 
critical point at 12T. The "critical point" is reinterpreted as a (noncritical) Kauzmann point in 
which the latent heat vanishes and the line is parallel to the T axis. The magnetization, the entropy 
Q ■ and the specific heat discontinuities at melting compare well with experiments. 
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The H — T phase diagram of a high T c superconductor is very complex due to the competition between thermal 
fluctuations (TF) and disorder. At low fields vortex solid melts into a liquid [1] due to TF. The discontinuity in 
magnetization [2-4] and in the specific heat [5,6] unambiguously demonstrate that the transition is a first order. 
Evidence is growing that the solid-glass transition in YBCO [7] is also a first order [4,8] . It was suggested [4] that the 
two first order phase transition lines H*(T) and H m (T) join at a multicritical point or that the melting line continues 
as a second order transition between the putative liquid II and a liquid [6] . However other experiments are interpreted 
using a concept of "unified first order transition line" [9]: only character of the transition evolves from the thermal 
q ■ fluctuation dominated to the disorder dominated one without multicritical points along the line. This single line was 
demonstrated in anisotropic materials BSCCO [10], LaSrCuO [11] and NdCeCuO [8] and was claimed recently in 
1 '! YBCO [12]. 

' Theoretically idealized models like the frustrated XY model [13,14] or a collection of interacting pointlike objects 
subject to both the pinning potential and the thermal bath Langevin force were simulated numerically [15,16]. Al- 
ternatively the melting line was located using the phenomenological Lindemann criterion on the solid side either 
in the framework of the "cage" model [9] or a more sophisticated approach based on the elasticity theory [17,18]. 
Although both approaches are very useful in more "fluctuating" superconductors like BSCCO, a problem arises with 
their application to YBCO close T c : vortices are far from being pointlike and even their cores significantly overlap. 
As a consequence the behavior of a dense vortex matter is different from that of a system of pointlike vortices and of 

' the XY model. 

In this Letter we quantitatively study the effects of the competition of thermal fluctuations and disorder in the 
framework of the anisotropic Ginzburg - Landau (GL) model appropriate precisely in the region of interest for YBCO. 
It was used successfully to describe the "clean" case [19,20]. 
The model without disorder is defined by free energy [1] : 
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where $o = hc/(2e), t = T/T c , A = (By, 0,0). The model provides a good description of thermal fluctuations as 
long as 1 — t — b << 1, where b = H/H C 2- An "applicability line" 1 — t — 6 << 0.2 is depicted on Fig.l. The model 
however is highly nontrivial even without disorder and within the lowest Landau level (LLL) approximation [19] in 
C$ ' which only the LLL mode is retained and the free energy simplifies (after rescaling) : 
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The simplified model has just one parameter - the (dimensionless) scaled temperature: ax = (t+b—1) <t/2/(Giir 2 b 2 t 2 ) 
with the Ginzburg number Gi = 32 (ttA 2 T c 7/($^))' ', j 2 =m c /m a i, the anisotropy parameter, A the magnetic penetra- 
tion depth and £ the coherence length. The (effective) LLL model is applicable in a surprisingly wide range of fields 
and temperatures determined by the condition that the relevant excitation energy e is much smaller than the gap 
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between Landau levels e c = 2heB/(cm a b)- Within the mean field approximation in the liquid e/e c — x is a solution of 
the "gap equation" [21]): 2bx + (1 - t - b) = 2irtVGib( v / 4x - 2 + J2n=aW 4x + 4n + 2 - y/4x + 4n - 2 - 1/^/n + x)]. 
The LLL dominance line in Fig.l represents a conservative condition e/e c = 1/20. On the right of the line, the HLL 
(higher Landau level) TF modes should be considered, while for b < (1 — 1)/13, HLL in solid becomes important [20]. 
Apart from the fields smaller than Hlll ~ 3fcG, experimentally observed the melting and the solid-glass lines are 
well within the range of applicability of the LLL approximation. 

The disorder is described within the GL approach via a random field W x with correlator W x W y — WS(x — y) [1], 

while TF are described by partition function Z = J^, ^ exp j — F + J x W x |iAz| 2 J /r j. Assuming small W, the free 
energy —TlnZ is calculated perturbatively to the second order: 

g = 9d + Jw x ^ x \ 2 )- ^^^[(|^| 2 |^| 2 )-(|^| 2 )(|^| 2 )]|, (3) 

where () and g c i denote the thermal average and effective free energy of the clean system. Averaging random field 
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and n{t) oc W/T. 



one obtains in the scaled units g = g c i — n{t)Agdi S , where Agdis = \ j 

Recently a quantitative theory [20] of TF in "clean" vortex liquids and solids based on the GL model of eq.(2) 
was developed and was successfully applied to the fully oxidized YBa^CuzOj for which the effects of disorder are 
negligible in whole currently accessible fields range [1,4]. The two loop result for solid in the clean case g so i — 

~a T /(2(3A) + 2.85 IotI^ 2 + 2.4/ot is sufficiently precise (on the 0.1% level) all the way to melting. We calculated the 
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one loop disorder correction: Ag so i = 2.14 \ar\ ■ An explicit expression for gu q (aT), obtained using the Borel-Pade 
resummation of the renormalized high temperature series is rather bulky and can be found in [20]. The disorder 
correction in liquid can be obtained by differentiating the "clean" partition function with respect to parameters: 
Agu q = (guq — 2aTg' liq )/3 — (g' liq ) 2 /2. These results enable us to find the location of the transition line and calculate 
discontinuities of various physical quantities across the transition line. 

It was noted [20] that in a clean system a homogeneous state exists as a metastable overcooled liquid state all the way 
down to zero temperature. This is of importance since interaction with disorder can convert the metastable state into 
a stable one. Indeed generally a homogeneous state gains more than a crystalline state from pinning, since it can easier 
adjust itself to the topography of the pinning centers. Since at large \clt\, Agu q oc is larger than Ag so i oc laT^ 2 , 
the transition line shifts to lower fields. The equation for the melting line is d(ar) = (guq— 9soi)/{Agu g — Ag so l) = n(t). 
The universal function d(aT), plotted in Fig. 2 turns out to be non-monotonic. Since n(t) is a monotonic function of 
t, one obtains the transition lines for various n in Fig.l by "sweeping" the Fig. 2. A peculiar feature of d(a-r) is that 
it has a local minimum at cit ~ —17.2 and a local maximum at ax ~ —12.1 (before crossing zero at ax ~ —9.5). 
Therefore between these two points there are three solutions to the melting line equation. As a result, starting from 
the zero field at T c , the transition field H(T) reaches a maximum at E ent (point CP in Fig.l) beyond which the curve 
sharply turns down (this feature was called "inverse melting" in [10]) and at E mag backwards. Then it reaches a 
minimum and continues as the solid (Bragg glass) - vortex glass line roughly parallel to the T axis. 

The temperature dependence of the disorder strength n(t), as of any parameter in the GL approach, should be 
derived from a microscopic theory assuming random chemical potential or fitted to experiment. We conjecture (and 
find to be consistent with the experiments) that the general dependence near T c is: n(t) = n(l — t) 2 /t (or W oc (1 — t) 2 ). 
The expression approaches the one used lower temperatures [1] with n = 77^ / (in 2 \J 2Gi^ 3 ) . The best fit for the low 
field part of the experimental melting line H m (T) of the optimally doped YBCO (data taken from [6], T c — 92.6, 
7 = 8.3) gives Gi = 2.010~ 4 , H c2 = 190 T, k = A/£ = 50 (consistent with other experiments [4,22]). This part is 
essentially independent of disorder. The upper part of the melting curve is very sensitive to disorder: both the length 
of the "finger" and its slope depend on n. The best fit is n = 0.12. This value is of the same order of magnitude as the 
one obtained phenomenologically using eq.(3.82) in [1]. We speculate that the low temperature part of the "unified" 
line corresponds to the solid - vortex glass transition H*(T) observed in numerous experiments [4,6-8,12] (squares in 
Fig.l [4]). A complicated shape of the "wiggling" line has been recently observed [12]. 

Magnetization and specific heat of both solid and liquid can be calculated from the above expressions for free 
energy. Magnetization of liquid along the melting line H m {T) is larger than that of solid. The magnetization jump is 
compared in Fig 3a with the SQUID experiments [5] in the range 80 — 90K (triangles) and of the torque experiments 
(stars [3] and circles [4] ) . One observes that the results of the torque experiments compare surprisingly well above 
83K while the SQUID data a bit lower than theoretical or the torque one. But those of [4] vanish abruptly below 
83K unlike the theory and are inconsistent with the specific heat experiments discussed below [6,22]. We predict 
that at lower temperatures (somewhat beyond the range investigated experimentally so far) magnetization reaches its 
maximum and changes sign at the point E mag (at which magnetization of liquid and solid are equal) . 
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Entropy jump calculated using the Clausius - Clapeyron relation is compared with an experimental one deduced 
from the spike of the specific heat [6] (triangles in Fig. 3b) and an indirect measurement from the magnetization 
jumps [4] (circles). At high temperatures the theoretical values are a bit lower than the experimental and both seem to 
approach a constant at T c . The theoretical entropy jump and the experimental one of [6] vanish at E ent (Fig.l) near 
75K . Experimentally such a point (called Kauzmann points [23]) was established in BSCCO as a point at which the 
"inverse melting" appears [10]. Below this temperature entropy of the liquid becomes smaller than that of the solid. 
Note that the equal magnetization point E mag is located at a slightly lower field than the equal entropy point E ent . 
The Kauzmann point observed at a lower temperature in YBCO in [8] is different from E ent since it is a minimum 
rather than a maximum of magnetic field. It is also located slightly outside the region of applicability of our solution. 
The point E ent is observed in [12] in which the universal line is continuous. 

In addition to the spike, the specific heat jump has also been observed along the melting line H m (T) [5,6,22]. 
Theoretically the jump does not vanish cither at E ent or E mag , but is rather flat in a wide temperature range. Our 
results are larger than experimental jumps of [6] by a factor of 1.4 to 2 (Fig. 3c). In many experiments there appears 
a segment of the second order phase transition continuing the first order melting line beyond a certain point. In [6] it 
was shown that at that point the specific heat profile shows " rounding" . We calculated the specific heat profile above 
the universal first order transition line. It exhibits a "rounding" feature similar to that displayed in [6] with no sign 
of the criticality. The height of the peak is roughly of the size of the specific heat jump. We therefore propose not to 
interpret this feature as an evidence for a second order transition above the first order line. 

As seen from the fitting the optimal doping YBCO, the disorder parameter n in the cases of interest is very small. 
However to address a more delicate question of appearance of the glass transition between liquid and glass claimed 
in some experiments [4,6] and simulations [15,14], we used a nonperturbative method, the replica trick, combined 
with the gaussian variational approximation [24]. Breaking of the replica symmetry signals appearance of the glass 
transition in static phenomena. Here we summarize the results, leaving details for a longer presentation. The most 
general hierarchical homogeneous (liquid) Ansatz [24] was considered and found that there is no replica symmetry 
breaking (RSB) solution. 

Now we compare our results with other theories starting with models based on the Lindemann criterion [9,17,18]. 
Near T c (where effects of disorder are small) the location of the first order transition line H(T) is qualitatively 
consistent with that found from the Lindemann criterion [1]. In the intermediate region around E ent our results for 
the melting line are completely different from all the phenomenological models [9,17,18]. In particular we do not have 
a critical point assumed in some of them (the two papers in [17] differ on the continuation of the line beyond the 
critical point in which two first order lines join together). At temperatures below 70i^ the elasticity theory approach 
based on the London approximation [17,18] (valid beyond the range of applicability of the GL approach) is expected 
to smoothly interpolate to the GL approach. The advantage of the present approach is that, in addition to location 
of the transition line, it enables the determination of discontinuities and the calculation of physical quantities away 
from the transition line. 

The comparison with numerical simulations can be made only qualitatively since the disordered GL model has 
not been simulated. In the XY simulations of [13] a single transition line is parallel to the T axis below certain 
temperature, while ours is not. In [14] there is, in addition, a slush - liquid transition, while Langevin simulation [15] 
finds the liquid - glass transition. Their transition criterion however is dynamical. Absence of RSB in liquid does not 
generally imply that this state, especially at low temperatures, does not exhibit "glassy properties" in dynamics. In 
fact overcoolcd liquids generally are "glassy". Therefore we propose to consider the glassy properties of the vortex 
state above the H*(T) line in the context of the "disordered overcooled" liquids. Recent simulations [16] demonstrate 
that glassy features in dynamics do not necessarily correspond to the conventional vortex glass scenario. Next we 
comment on the "vortex loops" scenario for YBCO. The region in which the loops are relevant according to [25] is 
below GiH C 2, too small to explain the experimental transitions (our fit for Gi = 2 10~ 4 is much smaller than that 
assumed in [25]). The LLL is inapplicable in this region (see Fig.l). In strongly anisotropic materials, a model of 
the Lawrence - Doniah type is more appropriate [1,26]. We performed a 2D GL calculation and found that in that 
case there is no "wiggle" of the transition line and speculate that it disappears at certain value of the anisotropy 
parameter. 
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FIG. 1. Theoretical first order transition lines for various degrees of disorder separated between a homogeneous and a 
crystalline phases of vortex matter. The best fit n — 0.12 line is compared with experimental melting line H m (T) ( [6]) and the 
solid - vortex glass transition line H* (T) ( [4]). Equal entropy and magnetization points are denoted by _E ent and E mag . The 
range of applicability is restricted by the two lines. Very small fields below Hlll are beyond this range. 



FIG. 2. Non-monotonic dependence of the disorder function d(ar) on the LLL scaled temperature ar 

FIG. 3. Discontinuities along the first order transition line in YBCO. Jumps of magnetization, entropy and specific heat are 
compared with experiments in a, b, c respectively. 
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